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ABSTRACT 


Observational evidence suggests that magnetospheric substorms may be 
associated with the formation of a pair of n utral points or lines in the 
geomagnetic plasma sheet, containing an X-type point (or line) and an 
0-type one. While magnetic merging theory has concentrated almost entirely 
on X-type neutral configurations (points, lines or sheets), here the role 
of 0-type configurations is examined, with special attention to three points: 
(1) How does the X-0 configuration extend in three dimensions? To this end, 
an analytical model of the configuration was derived, useful for visualizing 
the geometry and for numerical treatment of plasma flows in it; (2) What 
modifications are needed in the MHD condition E = v x B near the 0-type 
line, where it tends to make ^ grow without limit? By analyzing equations 
of motion ,for charged particles near an 0-type neutral line and their solu- 
tions in limiting cases, it was found that at a certain distance from the 
neutral line the mean particle motion became decoupled from that of magnetic 
field lines (which obey the MHD condition). The decoupling distance depended 
on initial conditions in momentum space, suggesting that the MHD approxima- 
tion which averages out such conditions may not suffice for describing 
plasma dynamics near the neutral line. Similar problems arise with merging 
flows near X-type neutral lines, and although the treatment there is more 
difficult and requires more approximations, it appears that the same qualita- 
tive conclusions apply there as well; and (3) IVhat is the role of 0-type 
neutral lines in particle acceleration? It was found that after inflowing 
particles are decoupled from the field line motion, they go over to a mode 
of runaway acceleration along the neutral line. This process is much more 
efficient along an 0-type line than along an X-type line and it is concluded 


1 


that if merging occurs at an X-0 pair, two particle populations may be 
expected — low energy particles accelerated adiabatically by earthward 
convection past the X-type line, dependent mainly on the total amount of 
flux which has been merged, and high-energy particles convected towards 
the 0-type line and undergoing there runaway acceleration. The second 
acceleration process depends critically on the rapidity of merging and 
is therefore expected to vary considerably from event to event. All 
this agrees with observations, and similar processes may also be 
important in solar flares, where a -'Y-type neutral point" has been 
proposed, which actually represents the limiting form of an X-0 
configuration. 
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CLASSIFICATION OF NEUTItAL POINTS 

The structure of magnetic field lines near a simple neutral point 
was first investigated by Dun gey [1953, 1963] who based his work on the 
Taylor expansion of the magnetic field in the vicinity of some given 
point P. This expansion begins x^ith 

® = ?0 

where B and VB^ are computed at the point P and the radius vector r 
is drawn from that point. If P is a neutral point then ^ vanishes and 
along any field lines passing through P the magnetic field must satisfy 

B = X, r (2) 

Comparison with (1) shows that must be a real eigenvalue of the 
dyadic and that ir is aligned with the corresponding eigenvector. 

Dungey noted that since may have either 1 or 3 real eigenvalues, 
there should exist two main classes of neutral points, having either 
1 or 3 field lines passing through them. He named such points 0-type 

and X-type neutral points, respectively, since their expected configura- 

* 

tions resembled somewhat the letters 0 and X (Figures 1). 

Dungey’ s analysis was repeated in a rigorous fashion by Fukao 
et al . [1975] who showed that the field line structure of an 0-type 
point consists in general of a set of spirals rather than of nested 
closed curves; however, if there exists a neutral line along which 
^ = 0, the field line structure indeed resembles Figure lb. Their 
analysis covered all possible configurations and included the special 
cases of degenerate eigenvalues and eigenvalues which vanish. 
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In most studies involving neutral points in space attention has 
been focused on X-type neutral points or neutral lines, or on their 
limiting form, neutral sheets. The main reason for this preference 
was that X-type points can always be expected to occur on the boundaries 
of different topological regimes of field lines emanating from two or 
more sources f Sweet , 1958] , A prime example of this Is the "open 
magnetosphere" model [Dun gey , 1961] where at least two X-type points are 
always expected on the boundary between open and closed field lines 
(Figure 2a). In contrast, no 0-type points are required in such a 
model. 

A further reason for the interest in X-type neutral points or 
lines has been that they form a key ingredient in theories of magnetic merging o 
field lines originating from different sources and embedded in a plasma 
in which the MHD approximation 

E = - V X B (3) 

is valid [e.g. Vasyliunas . 1975]. Because such a merging process can 
transfer magnetic flux from one topological regime to another and over- 
come the "freezing" of flux to the plasma [e.g. Stern , 1966] it has been 
widely believed that merging is responsible for solar flares and 
magnetospheric substorms, where magnetic energy appears to be rapidly 
converted and where charged particles are accelerated to high energies. 

Indeed, observations near the midplane of the geomagnetic tail have 
shown that while at geomagnetical ly quiet times the north-south magnetic 
component there is generally directed northward, as expected for 
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disturbed dipole field lines, during substorms it may briefly reverse 
to become southward [e.g. Hones , 1975, 1977; Nishida and Hones . 1974], 
suggesting the existence of X-type structures. 

O-'fYPE STRUCTURES IN THE TAIL 

Originally it was believed [e.g. Axford et al., 1965] that the 
neutral point involved in substorms was the nightside X-type point of 
the open raagnetospheric configuration of Figure 2a. However, more 
recent observations suggest that this is not the case — that while the 
raagnetospheric tail extends well beyond the moon's orbit (60 Rg) before 
becoming pinched off, substorra activity originates much closer to earth, 

at distances of the order of 10-15 R_. 

E 

Thus, if an X-type neutral point is involved in the mechanism of 
substorms, it arises in the interior of the tail, in a configuration 
resembling Figure 2b [Schindler . 1975; Hones, 1977]. As can be seen, 
this configuration also requires the existence of an 0-type neutral 
point or line. 

The purpose of this work is to use analytical models in order to 
examine three questions related to such 0-type neutral structures 
(points or lines) ; 

(1) What is the 3-dimensional form of the field in Figure 2b? 

(2) According to a commonly expressed view, when reconnection occurs 

in an X-type structure in the geomagnetic tail, a certain amount of magnetic 
flux originating in the northern lobe (in Figure 2b) flows through the 
structure, merges with an equal amount of flux from the southern lobe 
and then the combined field lines convect earthwards. As they con act, 


such field lines becomt> r.iort hi i olo-l ihe , reducing the stress in the 
magnetosphere and releasing magnetic energy. 

However, if an 0-typo line exists as in Figure 2b, an equal 
amount of flux flows into the "bubble" surrounding the 0-type line, 
converges towards the line and finally vanishes there. Does such 
behavior contradict MHD theory? 

(3) Finally, the role of an 0-type field configuration in the 
acceleration of charged particles is examined. As will be seen, there 
exist reasons to believe that particle acceleration can proceed much 
more efficiently inside the "bubble" than in the X-type region. 

In what follows, these three questions are addressed in the order 
given above. To simplify the discussion it will be assumed that we are 
dealing with 0-type lines similar to the one shoxm in Figure lb. A 
more general definition of an 0-type line which is sometimes used 
(e.g. Vasyliunas [1975], p. 304, penultimate parag.) allows a small 
component of B to persist along the "neutral line" Cexcept at isolated 

fs,* . 

"neutral points" where = 0): according to Fukao et al . , [1975 ], the 
field near that line would then have a small radial component, which 
would transform the nested closed field lines into spirals. Such 
fields are much more difficult to analyze and we presume that the 
components thus added are small and that the conclusions reached here 
are still qualitatively correct. 

A NALYTICAL MODEL FOR "BUBBLE IN TAIL" 

The only attempt to describe how the interior X and 0 lines of 
Figure 2b might be extended to the third dimension was made by Vasyliunas 
[1976, Figure 5 ] . Vasyliunas suggested that initially the 0-type line 


and the X-type line in the middle of the plasma sheet were topologically 
continuous. During a substorm, he proposed, both expanded until they 
reached the flanks of the tail, at which point the magnetic linkages 
became shuffled, the X^type line wound up connected to the dayside 
'toungey" point (Figure 2a) while the 0-type line linked up with the 
nightside neutral point. In what follows we shall investigate an 
analytical model of the initial configuration suggested by Vasyliunas, 
which is practically dictated by the over-all topology; the later 
developments mentioned above are much more speculative and are not 
relevant to the present sutdy. 

The analytical model is produced by taking a simple 2-dimensional 
model of the tail field 

B = (z/L) 6 X + z (4) 

(Bq and L are constant, the x direction is sunward and z = 0 is the 
midplane of the tail) and superimposing upon it a current ring in the 
z = 0 plane, flowing counterclockwise when viewed from z > 0, The aim 
here is pure electromagnetic modeling and no attempt is made to derive 
a plasma population which supports such a field in a self-consistent 
manner. 

Jackson [1962 ] has provided formulas for the magnetic field 
of a circular current filament of radius a, but these do not suit the 
present purpose, since b' becomes infinite as one approaches the fila- 
ment. What is required for a realistic model is a ring carrying a 
distributed current density varying smoothly enough so that when it is 
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added to the field (4), 15^ is reversed near tlie middle of the rliiR but 
j nowhere beeomes unrealistically large. 

Fortunately, Jackson also gives a simple approximation b which 
tends to b' at large distances and also near the origin and near the 
z~axis, and which fulfills all the above requirements. Tluit field 
given fin spherical coordinates r, h, ) by the vector potential 

■50 -a/9 

A = "f b_ r slnU fa^" + + 2ar sino) ‘ (5) 


leading to field components 
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C6a) 


(6b) 


Tlie above approximation represents a broad curret. t ring producing 
field lines which enclose the ring r = 2a in the z = 0 plane (although 
at large distances they approximate the field of a filament at r - a). 

In the z - 0 plane of the corresponding cylindrical system (i.e. the 
plane u = '■ /2) b^ Is southward (I.e. negative) wfienever r < 2a and 
northward (i.e. positive) whenever r 2a, but the magnitude j is 
small everywhere in that pl.are except near the origin (Table 1), so 
that when b is added to the field of equation (4), the reversal region 
is small, with a diameter of the order of a (while b reverses over a 
diameter 4a). It i.s surrounded by true X~type and 0— type neutral lines, 
since both b of (6) and B of (4) have only z components in the plane 
z ~ 0 and wherever these components cancel the total field vanishes. 

S 
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Before calculating and plotting tie combined field it is useful 
to tran^.fo^m everything into normalized units — i.e. scale all field 
intensities in units of Bq and all distances in units of the scale 
length L. Tlie configuration then depends on two dimensionless parameters — 
the normalized intensity of the ring current field 


C 
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(7a) 


and the no.-'inalized radius of the ring 

= a/L (7b) 

The range of these parameters is rather limited. The first para- 
meter must be negative and less than -0.5, otherwise no reversal 
of B occurs, and it is not Jikely to fall below -2, since then the 

southward B peaks at 3 B_, whici is about the limit of observed south- 
z u 

ward fields. Since L in the tail is of the order of 1 Rg, will 
approximately equal the diameter of the reversal region in earth radii , 
which can be taken to be around 1. 

Figures 3 show results for = -1.5, = 2. Specifically, they 

present projections of field lines onto planes of constant y/L (where 
y is the dawn-dusk coordinate), corresponding to y/L = 0, 0.25, 0.5, 0.55 and 
0.75. As can be seen, the neutral lines approach each other as jyj 
increases and they vanish at a value of y/L slightly exceeding 0,5. 
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The value oT B,, wa.^' alao computed; was found 

to he rather small, ;;u«ites Cliu: that to a rather good approximation 
the field lines do not depart far from planes of constant y. 

Topologically tlu model contains a distinct region, a "bubble" 
centered on the 0-type line, in which field lines are not connected to 
the main geomagnetic field: its shape resembles somewhat that of a 
football CAmeriean style). Since the detailed analysis of particle 
dynamics in this region is difficult, it will be handled in the next 
section for a somewhat simplified geometry, in which the "bubble" Is 
straightened out to a cylinder, 

PARTICLE MOTION NEAR AN LINE 

As a model for the magnetic field in the vicinity of an 0-type 
neutral line we consider a field 

B = va X Vp (8) 


where in cylindrical coordinates (^, z) the Euler potentials (u-, p) - 
are given by 


= Bq (\’g/2 -o f>?0 (9) 


and 


P = X (10) 

(in the context of the geomagnetic tail It should be noted that the 
directioa taken here as the z axis is aligned with the solar magneto- 
spheric y axis, so that the system of coordinates used here is not 
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aligned with the one used in the preceding section) 


For ii 


B . B. S' 


while for y< the magnitude of the field decreases In proportion to 

• 9 

f and vanishes along the z axis. 


For simplicity^ the calculation below asisumes that 1 is always less 


than the external region ? > 9q could have been included, but this 


would have only lengthened the calculation without introducing any new 


qualitative feature. 


Assume now that an electric field 


E = E„ z 

nx. U BX* 


exists in this region and that the hulk velocity of the plasma obeys 


the MHD condition 


E = - V X B 


Since E and B are orthogonal, v satisfies 


V E X B/B 


which is the well-known relation for the elect- ic drift velocity. In 


the given configuration this v is everywhere directed inwards, towards 


the neutral line, and its magnitude tends to infinity as ^ 0 and 


B 0. Thus the velocity of magnetic field lines described by (3) 


converges towards the 0-type neutral line and its magnitude tends to 


infinity just before those lines vanish into oblivion at 0. 
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If V Is viewed as a field line velocity then such a behavior 
creates no difficulties, because the motion of magnetic field lines is a 
mathematical abstraction and is not physicallv measurable. On the other 
hand, such behavior would not be appropriate for the bulk velocity of the 
ambient plasma, which is an observable material velocity. 

If the ambient plasma obeys the guiding center approximation at 
distances of the order or greater, than far away from the neutral 
line the velocity v given by (13) will indeed be its bulk velocity, and 
the plasma flow there will converge inward with a gradually increasing 
velocity. Near ? = 0, however, the guiding center approximation breaks 
down and one may expect that the MHD relation (3) is no longer valid. 

In what follows it will be shown in more detail how exactly this transi- 
tion occurs, how the motion of the plasma becomes decoupled from that of 
magnetic field lines and what the consequences may be. 

The approach adopted is a slight generalization of the one described 
by Stern [1975] . Under the conditions described above, the Hamiltonian 
for a particle of charge q and mass m is 


^ = h ^ ^ ^Pz - - q 2 (U) 

where (py , p^, P^) Bre canonical momenta conjugate to (?,'f , z) and 
where 

= = - 

Suppose first that = 0 (i.e., no electric field exists). Then 
p and p are both constants of the motion, with values which will be 
denoted by g and h, and the entire motion reduces to a motion in a one- 
dimensional potential V 
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H = p^/2m + V 


V = (l/2m) [g^/s^ + (h - qco^)^] 


The potential V depends only on u = 9 , is non-negative and 
(except for the singular case g = 0, which must be treated separately) 
tends towards infinity both as ? 0 and as ? . Tlie particle is 

thus trapped in a potential well and the bottom of this well is located 


at some value of ? , satisfying 


3V( ?^)Du =■ 0 (17) 

The particle's motion oscillates in that potential well around 
so that ? may be viewed as representing a "generalized guiding center" — 
and indeed, for low energy particles far away from the z axis, it does 
tend towards 9 of the guiding center. 

Now let the electric field be Included, so that H is given by (14). 
ITie canonical equations of motion then give 


Pz = ^ ^0 ■ *^0^ 


(18a) 


m z = qCEgCt - tg) - c f ] 


The complete motion cannot be solved, but in principle one could 
derive an approximate solution using the approach of Sonnerup [1971 ], 
as follows. The limit Eq =0, described by (15) and (16), (s one- 
dimensional and therefore formally soluablei one can thus (in principle) 
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Utn'ive for it an adiabatio invariant 


.1 - ^ Pp df - (2m) (W - V(? d? 

= J(W, g, p ) (lOa) 

o 

where W is the constant of energy- If is no longer zero but remains 
relatively small, J is still an approximate constant of the motion: 
inverting (19a) to the form 

p = f(J, W, g) (19b> 

z 

and substituting (18a) gives then 

q Eg (t - tg) = f(J, W, g) (20) 

If J and g are assumed to be conserved, the last relation describes 
(at least implicitly) the variation of the energy W with time- 

For particle motion near a neutral sheet, J was derived explicitely 
in terms of elliptic integrals by Sonnerup [1971] , who then incorporated 
the effects of a weak electric field in the manner described here. For 
the present calculation, the analytical form of J is too complicated 
and we shall contend ourselves witli a semi-qualitative treatment. 

Stressing the limiting properties of the motion very far from the line 
z = 0 and very close to it. 

We begin by substituting (18a) into (14). The Hamiltonian then 
still has formally the form (15), with the potential 

V( p, z, t) = (l/2m) tg^/p- + q^[Eg(t - tg) - r ^ f j - q Eg z (21) 
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However, z is now no Uniger a canonical coordinate, since is gone: 
instead. It is parameter evolving according to (18) (the canonical 
relations for 9 and remain the same as before). For any fixed z 
the dependence of V on 9 is as before a potential well bottoming at 

f) 2 

\ , representing a "generalized guiding center". If u ~ 9^^ then by 

07) 

3V/3U = K(ii) - - a“/u^ + 2 [E^Ct - t^) - cu ] = 0 (22) 

The difference "ow is that due to the dependence on t, 9^^ 
gradually shifts, i.e. 3 ^ differentiating (20) 


from which 


2 ?j^(3 9j^/at) (3K/3u) -f K/3t = 0 



2 ? - (2g^/q^ 


(23) 


(24) 


This result is easily interpreted. Far from tlie z axis, where 
is large, the second term in tl:e denominator is small and may be 
neglect(>d. By equations (8) - (10), the first term there is simply 
-B, so the flow in this limit is directed radially inwards with the 
electric drift velocity Eq/B: the mean motion in this region is thus 
in accord with (3) and (13). 

However, at some critical distance 



= (-g^/q^c) 


1/6 


( 25 ) 


(note Chat c is negative) the tA-ro demoninator terms become equal, while 
beyond this distance the second term dominates. The motion is then 
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still inwards (since both terms are negative) , but the rate at which 
P decreases slows down to a negligible value. 

Ifliat happens to the particles? If Pj^ is viewed as a "generalized 
guiding center",, the motion can be separated into a mean part and an 
oscillating one, and in particular it is possible to write 

p' - + ( 26 ) 

For studying the average properties of the motion in the z direction, 
the oscillating part in (18b) may be neglected, giving after differentia- 
tion 

" - 2 <1 fl (2 ?i/at:) (27) 

For large values of ^ by (24), the two terms tend to cancel and 
to the lowest approximation no appreciable acceleration takes place. A 
more accurate calculation would show that the particle actually loses 
some energy, because in the outer regions the magnetic moment is 
approximately conserved (this conservation may be viewed as the limiting 
case of (19a)). As decreases and the particle moves Inwards, it 
enters regions of progressively weaker B and if |a is preserved, its 
energy W also decreases: however, the electric drift velocity, which 
dominates its motion in these regions, does not depend on W, so that 
the particle's inward progress continues undiminished. 

Once has decreased to the order of or less, however, the 

second term in (27) loses importance and the particle is accelerated almost 
freely along the z axis. These two regimes are completely analogous to 
corresponding regions derived by Sonnerup for the vicinity of a neutral 
sheet. 
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The apparent conflict with MHD is thus resolved as follows. 

Magnetic field lines in the 0-type configuration flow radially Inwards 
with a steadily increasing velocity given by (13), and they vanish at 
z = 0 like the mythical snake which swallows its own tail. This 
velocity however is not physically observable since B In that region 
(though not at great distances, where equations (8) - (10) no longer 
hold) is time Independent. 

Plasma particles which start out by sharing the inward flow of 
magnetic field lines decouple from the magnetic field flow at distances 
of the order of ^^d depends, for each particle, on the Initial 

conditions (through the parameter g). Once the decoupling has occurred 
the particles are removed from the scene by being freely accelerated by 
E along the middle of the 0-type configuration, which acts as an 
accelerator tube. 

The next section examines the extent to which this kind of 
"decoupling" occurs in the neighborhood of an X-type neutral line. 

X- TYPE NEUTRAL LINES 

One of the major problems in merging theory has been that if the 
MHD condition is assumed to hold in the merging plasma 


E - - v X B (3) 

then some modification must be added to it before it can be applied near 
B = 0. Such a modification has been frequently introduced by regarding 


(3) as an approximation to the generalized Ohm's law equation [e.g. 
Vasyliunas , 1975, eq. 1"; 1976, eq* 2; Rich et al., 1976, appendix] 
and retaining some other terms from this equation near the neutral 
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point — for Instance, the resistive term (with normal or anomalous 
resistivity) or inertia terms. 

In the preceding section, in connection with merging at an 0-type 
neutral line, a different approach was adopted — namely, (3) was 
regarded as an approximation to the mean motion of a charged particle 
in the plasma. It was then shown that near an 0-type neutral line, 

(3) breaks down and a different regime of motion takes over, one in 
which particles are steadily accelerated. This breakdown cannot be 
easily accommodated by MHD theory, because it occurs at different 
distances for particles which differ by the constant of motion p^ , 
and such particles are not distinguished in MHD theory. 

It is natural to ask whether such a calculation can be extended 
to X-type neutral lines. Unfortunately, a new difficulty is encountered 
here; while the motion near an 0-type point in the absence of an 
electric field can be reduced to one-dimensional form (equations 15-16), 
the corresponding form for an X-type line remains 2-dimensional 
£Stern , 1975] and cannot be integrated analytically. The best one can 
do by analytic methods is to trace the motion of the guiding field line — 
rather than that of the guiding center — as shoiai below. 

If B is given by Euler potentials as in (8) and cartesian coordinates 
(x, y, z) are used, with the z axis along the neutral line, one finds 
here 

a = M^x^ - N^y^ (28a) 

= z C28b) 


18 


The field lines (Figure 4) are now a set of hyperbolas sharing 
the asymptotes 


MX = + Ny (29) 

which correspond to a = 0, and in the absence of electric fields the 
Hamiltonian is 

H « (l/2m) (p^ + p^) +V (30) 

X y 

V - (l/2m) [p^ - q(M^x^ - N^y^) f (31) 

Tlve "potential" V is non-negative, and since p is a constant of the 
motion, V reaches its absolute minimum value zero along the "guiding 
field line" 

a == p^/q = a.Q (32) 

At infinity (x ^ or y — oo) V also becomes infinite, while at the 

2 

origin it reaches a local maximum p^/2ra. Thus particles for which the 

(x, y) component W of the kinetic energy is sufficiently low, or 
xy 

which are sufficiently distant from the origin, will be confined to a 
narrow valley centered on (X = ttQ, and by analogy with the preceding 
section one can state that the particle's generalized guiding center 
is located somewhere along the bottom of this valley: it is not 
possible to give the position more precisely, since V is now two- 
dimensional. 
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Actually^ there exist two hyperbolic valleys along which V vanishes, 
located symmetrically on opposite sides of the origin. When the particle 
l^s far i'rora the origin, that valley In which the particle Is not located 
may be ignored. On the other hand. If the particle comes close to the 
origin and if W is sufficiently large to enable it to cross the pass 
at the origin (Figure 5), then can no longer be regarded as the 
generalized guiding line — instead, the particle now oscillates 
irregularly in the region surrounding the origin [Rus abridge . 1971, 1977], 
Let now an electric field z be added as before. In analogy with 
with (15) and (21) we now have 

H * (l/2ra) (p^ + p^) + V(x, y, z, t) (33) 

V = (q^/2m) [EQ(t - t^) - + N^y^]^ - 9 z (34) 


where z is to be regarded as a parameter satisfying the analog of (18b) 

m z = q [EQ(t - tp) - + N^y^ ] (35) 

For any constant z, the bottom of the potential "valley” is along 
the hyperbola 


mV - nV = EgCt - tg) (36) 

This can be regarded as the equation of the guiding field line 
a = which shifts in time towards increasingly larger values of a 


% - “ V 


(37) 


ORIGINAL PAGE IS 
GE POOR QUALITY 


Tills again reflects the electric drift velocity of the guiding 

center, given by (13). This velocity is known to be a valid field line 

velocity [e.g. Stern , 1966]; even if the position of the guiding center 

along the field line is not known, if each particle along the line moves 

with velocity (13), the same field line will continue to contain those 

particles at all subsequent times. If the energy V of the (x, y) 

xy 

motion is low, the particle will oscillate around its guiding field 
line and its a at any instant can be expressed as 

Substituting the last 2 equations and (26) in (33) gives 


(i) 


aver 


= 0 


(39a) 


and hence, in analogy with the result of the preceding section. 


( 2 ) 


aver 



(39b) 


In other words (as could be expected) the particle moves, on the 
average, with its electric drift velocity (13) and is not accelerated 
along z. In addition it may also move along its guiding field line 
and oscillate across it, and its energy changes slowly in accordance 
with the conservation of jo., but these details are not resolved by the 
present approximation. 

However, if the particle manages to approach the origin sufficiently 
close so that it crosses over from one valley to the other, all the 
preceding breaks down. Its average x and y now oscillate in an 
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irregular manner (investigated numerically by Russbrtdge ([1971, 1977]): 
it is not clear whether the average of (M x - try ) is zero, but on 
the other hand, as long as the particle stays near the origin, it 
remains bounded and it may be argued that when the average time 
derivative of (35) is computed, its contribution may be neglected. 

If this is the case then in the central region the particle again under- 
goes, on the average, a steady acceleration (as in the preceding case), 
satisfying 

^'^^aver = ^ 

Thus there exists a strong suggestion that particle motion near 
an X-type neutral line also contains two limiting regimes — an 
adiabatic drift regime far away from the line, in the region where the 
MHD condition (3) holds, and a runaway acceleration regime near the 
neutral line itself, where (3) breaks down. The existence of such 
regimes can be shown explicitly for one particular group of particles, 
namely those with initial conditions 

X = 0 p^ = 0 (41) 

From Hamilton’s equations, if (41) holds at one time, it holds 
for all other times as well, i.e. the particles are confined to the 
y axis. The fields sensed by them and the equations describinE their 
motion are then the same as those existing in Sonnerup's null sheet 
geometry fSonnerup , 1971 1 . where such limiting regimes can be derived 
analytically from adiabatic Invariance. 
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'llaere exists, however, one important difference between this case 
and the 0-type line. With the 0-typc configuration, every one of the 
particles swept inwards in the adiabatic mode ultimately passes into 
the second mode and becomes accelerated, unless E has meanwhile changed. 

rw 

With the X-type line the present treatment does not indicate what 
fraction of the particles attains the second mode and becomes accelerated, 
since all particles guided by the same field line arc handled together. 
Each such field line passes through the origin at t = t^ and at that 
time some of the particles associated with it will have passed to the 
second mode. Tlie point to note is that the fraction of such particles 
cannot be large, because if (13) is regarded as the bulk velocity in 
the first region, most trajectories of this motion will skirt the 
vicinity of the origin by a wide margin and only those among them 
which are headed almost directly towards the origin reach the region 
of acceleration (individual particle trajectories differ from the 
trajectories of the bulk motion by the addition of a component v„ 
parallel to magnetic field lines, but this does not affect the con- 
clusions reached here) . 

Furthermore, the duration of sustained acceleration in the vicinity 
of an X-type neutral line may be brief for most particles undergoing 
it, because the region in which it occurs has 4 "valleys" leading 
away from it. Numerical tracking of trajectories [Russbridge , 1971, 

1977] suggests that sooner or later particles enter these valleys — 
and although adiabatic mirroring will generally drive them back out 
again, their electric drift may meanwhile move them away from f le region 
in which (4Q) holds . 
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n>e relevant conclusion from all this Is that if during substorms 
two magnetic neutral lines of different types are created as shown in 
Figure 2b, close enough to each other to share the same electric field, 
then far more particles will undergo runaway acceleratio:i near the 
0~type line, which draws them in and keeps them well contained, than 
near the X-type line. 

S UBSTORltS AND FLARES 

lihat happens during a magnetic substorm? Observations seen to 
indicate fFairfield and Ness , 1970; Caan et al . , 1973] that stretched 
tail lines snap back to a more dipolar shape, decreasing the amount of 
stretched magnetic flux threading the cross-section of the near-earth 
tail, llie exact mechanism and cause of this leap are still a matter 
of controversy, but it is widely believed that it is initiated by the 
formation of an interior X-type neutral point or line as shown in 
Figure 2b, a view which will be adopted in what follows. 

The decrease in the stretched magnetic flux induces an e.m.f. in 
the circuit consisting of the plasma sheet and the magnetopause 
f Stern , 1977, sect. 6b] and this inductive e.m.f. is responsible for 
the enhanced electric field E generally assumed to exist along the 
neutral lines from da\«i to dusk. It should be stressed that E is a 
global phenomenon, due to the decreasing magnetic flux in the tall, and 
not a local consequence of conditions in the merging region: its energy 
is derived from the magnetic energy of the decreasing tail flux and it 
is available in its turn to accelerate and energize tall particles. 

Two simultaneous modes of acceleration may be expected to result 
from E. First, there would exist an enhanced convection of particles 

r<* 


streajning earthwards abound the X-llne: most of them will stay out of 
the immediate vicinity of the X-line, although their guiding field lines 
will pass through it. The energization of these particles will be 
moderate and it can be estimated from adiabatic conservation laws as 
being of the order of the ratio between the final and initial field 
intensities sensed by the convected particles. This process may be 
the one responsible for substorm particles in the ring current and the 
aurora, although the properties of such particles may tmdergo further 
modification in the inner magnetosphere by local electric fields with 
E„ i 0. 

Secondly, there exists a high energy component 0.5 Mev) detected 
in a number of substorms fRoelof et al .. 1976, Keath et al . . 1976; 

Hones et al . . 1976; Baker and Stone . 1977] and even observed outside 
the magnetosphere. It is proposed here that such particles are accelerated 
along the 0-type line of a configuration such as the one of Figure 2b, 

This line probably occupies a '’bubble" extending only 
part of the way across the tail and particles emerging from it become 
attached to magnetospheric field lines, which is the mode of motion 
in which the energetic particles are generally observed. 

In this connection It should be noted that there exist many 
similarities between substorms and ilar flares [e.g. De Felter , 1975] 
and that the radiation from flares suggests that there, too, the 
accelerated particles can be divided into two populations occupying 
different energy ranges fBai and Ramaty , 1976]. Tlie conventional 
explanation of this is based on a two-stage acceleration process — 
initially particles are accelerated to the lower level and a certain 
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fraction of them are accelerated again. If the present model is 
applicable, however, i.e. if flares occur at X-0 pairs of neutral lines 
within bundles of stretched flux emanating from sunspot regions <— these 
could be two distinct populations, accelerated in distant ways. In 
this connection it may be pointed out that the ”Y type neutral line" 
prc.posed by S tor rock [1973 ] closely resembles the limiting configuration 
of an X-0 pair when both of its components lie very close together. 
CONCLUSION 

The theoretical picture of magnetic merging and particle accelera- 
tion in substorms presented above is incomplete in many details. Not 
only can particle trajectories near neutral lines be handled analytically 
only in some limiting eases and in an idealized geometry, there also 
exists only a poor understanding of what precedes the substorra — of how 
the magnetosphere manages tc establish a configuration of plasmas and 
magnetic fields which can rapidly relax and release energy as soon as a 
pair of neutral points is established, but not before. 

We also lack information about the duration of the primary energy 
release in substorms (or flares), which is related to the rate at which 
magnetic flux jumps back. Stern [1977, sect. 6b] estimated that if the 
change associated with a typical substorm occurred over 30 minutes, the 
Induced e.m.f. would be 40,000 volts: if it takes only one minute, 
oi ir 10^ volts are generated and the energetic particles observed in the 
0,5 Mev range are readily produced along the 0-type line. 

In practice, the rate at which the flux changes may well depend on 
circumstances and vary from one event to the next one. If the present 
picture is correct, such variations will affect only slightly the energy 


transfer to plasma convected through tlie X-type line, for unless the 
depth of penetration changes for the convceted plasma, the same energy 
is provided to the same number of particles — only the time required 
tor the process differs. On the other hand, this variation will make 
a great difference in acceleration effects along the 0-type line. 

Substorms — and perhaps also solar flares — which lack a high energy 
component might well differ from those which possess it only in the 
rate at which they evolve, which might be slightly slower. 

To suiimiariae: 

(1) A model has been constructed for the magnetic field configuration 
of a combination of X-type and 0-type neutral lines in the middle of the 
geoniagentic tail, llie main feature is a "bubble” centered on the 0-type 
line, the field lines of which are not connected to the main geomagnetic 
field. 

(2) An analytical model of the motion of charged particles in the 
vicinity of an 0-type neutral line, with an electric field aligned 
parallel to it, has been constructed and analyaed. Using the concept 

of a generalized guiding center, it is shown that on the average particles 
move steadily inwards towards the 0-type line. At large distances their 
motion obeys the MHD relation, but as they pass within some critical 
radius the liHD approximation breaks do™ and particles tend to become 
freely accelerated along the 0-type line. 

Although the motion of charged particles cannot be analyzed as 
well in the neighborhood of an X-type neutral line, it seems that similar 
regimes exist there as well. There is a difference, however, in that 
only very few particles can here reach the region of free acceleration, 
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an*', ■'nfinement in this region is also less thorough than It is near 
an 0-type line. Tlje conclusion is Chat if particles are accelerated 
across the. plasma sheet along a neutral line, this is far more likely 
to be the 0-type line than the X-type line. It is also concluded that 
since nonadiabatic behavior is responsible for the breakdown of the 
MHD approximation near an 0-type neutral line (and probably near an 
X-type line as well), treatment of merging of collision- free plasma 
near such lines should not be based on the MHD approximation, even if 
eq. (3) is supplemented by resistive or inertial terms. 

(3) If a substorm is powered by the inductive e.m.f. generated 
when the tail’s magnetic field returns to a more dipole-like configura- 
tion by flowing through a pair of X-type and 0-type neutral lines 
(Figure 2b) , then two distinct mechanisms of pratlcle acceleration may 
take place. On the earthward side of the X-type line particles are 
moderately energized by convection: their final energy depends mainly 
on the depth to which they penetrate into the inner magnetosphere, while 
the total energy given to them depends mainly on the amount of reconnected 
flux. This process is not critically affected by the rate of merging. 

In addition, however, particles are also convected towards the 
0-type line and are accelerated there. Such particles can reach high 
energies, provided the rate of merging is sufficiently rapid. It is 
conjectured that while the first mechanism injects particles into the 
ring current and aurora, the second one is responsible for the fluxes 
of high energy particles recently observed during substorms by Imp 
spacecraft. Similar processes, occurring on stretched field lines above 
sunspots, may explain the existence of two populations of energetic particles 
in solar flares, commonly ascribed to a two-stage acceleration process. 
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FIGURE CAPTIONS 


Figure la - 


Figure lb — 


Figure 2a ~- 


Figure 2b — 


Figure 3a 


Figure 3b - — 


Figure 3c — 


Figure 3d — 


Structure of magnetic field lines near an X-type neutral 
point. 

Structure of magnetic field lines near an 0- type neutral 
point. This drawing illustrates a special case, corresponding 
to an 0-t3rpe neutral line along which B = 0. 

(fV 

Schematic view of the "open magnetosphere" proposed by 
Dungey. 

Magnetic merging at a pair of neutral points formed in the 
Interior of the tail, as proposed by some theories of 
magnetic substorras. One of these points has an X-type 
configuration, the other is of the 0-type variety. 

Magnetic field lines in the plane y = 0 of a three- 
dimensional model of the X-0 configuration, as described 
in the text. Distances are scaled in units of a, the 
characteristic dimension of the current loop, and the 
field structure for z < 0 is a mirror image of the one 
portrayed here. 

Similar lines for y = 0.25 (a/L). The small component 
Is ignored in this graph and in the three that follow 
it, so that the lines dravm actually describe two- 
dimensional projections of the magnetic field. 

Similar to 3b, but for y = 0.5 (a/L). The two neutral 
points are now rather close to each other. 

Similar to 3b, but for y = 0.55 (a/L). The lowest field 
intensity (at x = z = 0) is down to 0.0209 B^, but no 
neutral points remain anywhere. 


Figure 


Figure 


Figure 
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3e — Similar to 3b, but for y = 0.75 (a/L) . Tlie shape of the 
field lines is only moderately affected by the proximity 
of the current ring, but the field intensity at the 
origin is depressed to 0.3181 

4 — Schematic view of field lines in the x-y plane near a 

neutral line. Each line here can also be viewed as a 
line of constant potential , with V defined by 
equation (31) . The two dashed lines are the bottoms 
of two "valleys" along which V ■= 0, for some particular 
value of p : the choice of a different value v;ill shift 

the valleys to a different pair of equipotential lines. 

5 — The variation of the potential V along the line y * 0 

2 

in Figure 4. Here V is given in units of P^/2m (its 

2 2/2 

value at the origin) while x is in units of(p 2 /qM ) 
Particles in this potential may be trapped on one side 
of the origin or may cross from side to side: these 
modes correspond to two modes of motion distinguished 
in the work of Sonnerup (1971] and also represent two 
regimes of motion expected to exist in the vicinity of 
the X-type neutral line. 



TABLE CAPTIOM 

Table 1 — ITie values of B/Bq for the inagitetic field of equations (6), 

at various values of the radial distance r/a in the plane 
2=0, The scaling of distances is the same as in 
Figures 3; the field intensities along z = 0 in Figure 3a 
equal 1 - I.SCB/Bq), if one identifies x/a there with r/a 
in the table. Tlie value B/Bq = -1/256 at r/a = 3 is the 
most negative one encountered. 


TABLE 1 


r/a 

B/bo 

r/a 

B/b^ 

0 

2 

0.8 

0.114 

0.1 

1.30 

1.0 

0.0625 

0.2 

0.87 

1.5 

0.0128 

0.3 

0.60 

2.0 

0 

0.4 

0.42 

2.5 

-0.0033 

0.5 

0.30 

3.0 

-0.0039 
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